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MAGNETIZED NEUTRON STAR ATMOSPHERES: BEYOND THE COLD PLASMA APPROXIMATION 
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ABSTRACT 

All the neutron star (NS) atmosphere models published so far have been calculated in the "cold plasma 
approximation", which neglects the relativistic effects in the radiative processes, such as cyclotron emis- 
sion/absorption at harmonics of cyclotron frequency. Here we present new NS atmosphere models which 
include such effects. We calculate a set of models for effective temperatures T e s = 1-3 MK and magnetic 
fields B ~ 10 10 -10 11 G, typical for the so-called central compact objects (CCOs) in supernova remnants, 
for which the electron cyclotron energy E c , e and its first harmonics are in the observable soft X-ray range. 
Although the relativistic parameters, such as kT c s/m e c 2 and E c ^ e /m e c 2 , are very small for CCOs, the rela- 
tivistic effects substantially change the emergent spectra at the cyclotron resonances, E s» sE c>e (s — 1,2,...). 
Although the cyclotron absorption features can form in a cold plasma due to the quantum oscillations of the 
free-free opacity, the shape and depth of these features change substantially if the relativistic effects are in- 
cluded. In particular, the features acquire deep Doppler cores, in which the angular distribution of the emergent 
intensity is quite different from that in the cold plasma approximation. The relative contributions of the Doppler 
cores to the equivalent widths of the features grow with increasing the quantization parameter b e g = E C}e /kT e s 
and harmonic number s. The total equivalent widths of the features can reach ~ 150-250 eV; they increase 
with growing b c s and are smaller for higher harmonics. 



Subject headings: radiative transfer — stars: neutron — stars: magnetic fields — pulsars: 
(IE 1207.4-5209, PSR J1210-5226, PSR J1852+0040, PSR J0821-4300) 
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1. INTRODUCTION 

Thanks to the high sensitivity and spectral resolution of the 
X-ray observatories Chandra and XMM-Newton, absorption 
features in thermal spectra of the so-called X-ray dim isolated 
neutron stars (XDINSs 5 ) and central compact objects (CCOs) 
in supernova remnants (SNRs) have been discovered. The X- 
ray emission of these NSs is pulsed, and the measured period 
derivatives are consistent with strong (a few times 10 13 G) 
magnetic fields in XDINSs (Haberl 2007; Kaplan 2008) and 
relatively weak (< 10 11 G) magnetic fields in CCOs (Halpern 
& Gotthelf 2009). Therefore, the observed absorption fea- 
tures might be associated with ion cyclotron lines in XDINSs 
and electron cyclotron lines in CCOs. Correct identification 
of the observed features can provide important information 
on the chemical composition, magnetic field and gravitational 
redshift at the NS surfaces. 

The most famous isolated NS with absorption features is 
IE 1207.4-5209 (hereafter IE 1207), which is the CCO of the 
PKS 1209-51/52 SNR, with properties similar to other CCOs 
(Pavlov et al. 2002, 2004; de Luca 2008; Halpern & Gotthelf 
2009). Two absorption features in the spectrum of this iso- 
lated NS, centered at about 0.7 and 1.4 keV, were discovered 
by Sanwal et al. (2002) in Chandra observations, while later 
observations with XMM-Newton suggested the presence of 
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two more features, at 2.1 and 2.8 keV (Bignami et al. 2003). 
The period P = 0.424 s was found by Zavlin et al. (2000) 
from X-ray pulsations. IE 1207 shows a thermal-like X-ray 
spectrum with a black-body (BB) temperature ss 3 MK. The 
effective temperature T e g may be substantially lower (sa 1 
MK) if the spectrum is formed in a hydrogen atmosphere (Za- 
vlin et al. 1998). 

It has been suggested that the harmonically spaced absorp- 
tion features in IE 1207 are the electron cyclotron line and its 
harmonics (Bignami et al. 2003). At first glance, one should 
not expect spectral features at harmonics of the cyclotron en- 
ergy because in a classical (non-quantum) plasma the harmon- 
ics are due to relativistic effects (e.g., the ratio of emissivi- 
ties and opacities in consecutive harmonics is ~ kT/m e c 2 
at kT > E c<e ), which are very small in the relatively cold 
CCO atmospheres. However, Suleimanov et al. 2010 (here- 
after SPW10) have demonstrated that, even in a cold plasma, 
harmonically spaced absorption lines can form due to quan- 
tum oscillations in the dependence of the free-free absorption 
on photon energy in a magnetic field (Pavlov & Panov 1976; 
see also a recent discussion of this process in Potekhin 2010). 
The peak energies of these oscillations coincide with the cy- 
clotron energy and its harmonics, E = sE c _ e = sheB/m e c 
(s = 1,2,.. .). For IE 1207, they correspond to the mag- 
netic field Aine = 6 x 10 10 (^ e /0.7keV)(l + z) G (z is 
the gravitational redshift in the NS surface layers), consistent 
with B sc \ = 1 x 10 11 G, which corresponds to one of the 
two solutions obtained by Halpern & Gotthelf (201 1) for the 
pulsar's period derivative 6 . 

The model atmospheres were calculated by SPW10 in the 
so-called cold plasma approximation (Ginzburg 1970), which 
implies, in particular, that effects associated with thermal mo- 
tion of plasma particles can be neglected. All the previous 

6 This solution had been reported by Pavlov & Luna 2009 (see footnote 
#15 of that work). 
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magnetized NS atmosphere models, starting from the first 
work by Shibanov et al. (1992), have been calculated us- 
ing this approximation. Without allowance for the above- 
mentioned quantum oscillations, there is only one resonance 
(the fundamental cyclotron resonance at E = E ce ) in the 
electron component of the opacity of a fully ionized "cold" 
plasma, and this resonance occurs only in one polarization 
mode, called the extraordinary mode (X-mode), while the en- 
ergy dependence of the opacity of the ordinary polarization 
mode (O-mode) remains smooth (see Pavlov et al. 1980a; 
hereafter PMS80). The thermal motion leads to the cyclotron 
resonance (s = 1) in the O-mode (Pavlov et al. 1979), to emis- 
sion/absorbtion of the cyclotron harmonics (s > 1) in both 
polarization modes, and to the Doppler broadening of the cy- 
clotron resonances (e.g., PMS80); the contribution of these 
processes is determined by the ratio kT/m e c 2 . In addition to 
the thermal effects, quantum relativistic effects (whose con- 
tribution is determined by the ratio E c e /m e c 2 , can be impor- 
tant when E Cte > kT (Pavlov et al. 1980b; hereafter PSY80). 
These effects have been ignored in the NS atmosphere mod- 
eling because the parameter £ e = m&x(kT,E Cye )/m e c 2 is 
small (~ 10~ 3 ) in atmospheres of isolated NSs. However, £ e 
is not the only parameter that determines the strengths of the 
cyclotron resonances. For instance, the depth of the cyclotron 
absorption features depends on the ratio of the peak opacity 
to the opacity in the neighboring continuum, which may be 
^> 1 even for a very small £ e . Therefore, to properly interpret 
the spectral features observed in CCOs, one should use NS 
atmosphere models calculated with account for the thermal 
motion of electrons and quantization of the electron motion 
in a magnetic field. 

Here we present first calculations of realistic magnetized 
NS atmosphere models with these effects taken into account. 
We show that the fundamental resonance in the O-mode and at 
least two harmonics in both polarization modes are significant 
in the emergent spectra at the effective temperatures and mag- 
netic fields typical for CCOs. In particular, the thermal motion 
effects lead to formations of relatively deep Doppler cores at 
the centers of broad absorption features that arise due to quan- 
tum resonances in magnetic free-free absorption. These cal- 
culations demonstrate that the thermal motion and quantiza- 
tion must be taken into account for a quantitative description 
of the observed features in the spectra of CCOs. 

2. OPACITIES IN ELECTRON CYCLOTRON LINE AND HARMONICS 

Our consideration of the thermal cyclotron opacities is 
based on the works of PMS80 and PSY80, where simple ap- 
proximate expressions for the cyclotron opacities were de- 
rived. 

In PMS80, classic (non-quantum) relativistic opacities are 
considered with electron-ion collisions taken into account. 
According to that work, the shape of the cyclotron resonances 
can be described by the Voigt profile, which is the convolu- 
tion of the Doppler and Lorentz profiles. For ^> T c + T r 
(where T c , T r and To are the collisional, radiative, and 
Doppler widths of the cyclotron resonances), the Voigt pro- 
file consists of a Doppler core, independent of T c and r r , and 
Lorentz wings, which are virtually unaffected by the thermal 
motion. 

In the atmospheres considered in our work, the quantum ef- 
fects are significant because E c c > kT c s. The cyclotron 
absorption in the relativistic quantum case has been studied 
by PSY80 for a collisionless plasma. Since the free-free tran- 
sitions, caused by electron-ion collisions, are not included in 



that work (as well as in later works on thermal cyclotron pro- 
cesses in quantizing magnetic fields), the results of PSY80 
are applicable only in Doppler cores of cyclotron harmonics. 
However, similar to PMS80, one can use the cold plasma ap- 
proximation outside the Doppler cores (i.e., in the Lorentz 
wings and the continuum), because the relativistic correc- 
tions are negligible there (at kT e <C m e c 2 ). Moreover, since 
the thermal cyclotron opacity in the Doppler cores is much 
larger than the free-free opacity (except for the fundamental 
resonance in X-mode), and the thermal cyclotron opacity de- 
creases exponentially in the wings of the resonances, one can 
calculate the total opacity as the sum of these two (see Fig- 
ure 1 in PMS80). Therefore, we use the following approach 
in our calculations. For each used photon energy E, we cal- 
culate separately the opacities in the cold plasma approxima- 
tion, with the quantum oscillations in the magnetic free-free 
opacities taken into account (see details in Suleimanov et al. 
2009 and SPW10), and the cyclotron opacities (with quantum 
effects taken into account) in the collisionless approximation 
(PSY80). Then, we simply add the collisionless opacities at 
all the harmonics (except for the fundamental resonance in X- 
mode) to the magnetic free-free opacities in the corresponding 
modes. For the fundamental resonance in X-mode, we com- 
pare the two opacities and take the larger one as the actual 
opacity at this energy. The used approach only slightly dis- 
torts the frequency dependence of the absorption coefficients 
in the transition regions between the cores and the wings, but 
this should not significantly affect the emergent spectrum be- 
cause this transition occurs at different distances from the cy- 
clotron resonances at different depths. 

We perform all calculations in the local thermodynamic 
equilibrium (LTE) approximation. This means that the cy- 
clotron absorption is considered as a "true absorption", with- 
out scattering, and the cyclotron absorption does not con- 
tribute to the scattering part of the source function in the radia- 
tive transfer equations, which take into account the coherent 
electron scattering in the continuum only. 

The approximate equations that we use to compute the cy- 
clotron opacities are presented below (see PSY80 for details). 
They are not sufficiently accurate at angles 9 between the 
magnetic field and the photon wave vector close to 90°, i.e., 
at | cos 9\ < f3, where (3 is the ratio of the thermal velocity vt 
to the speed of light: 



c \m e c 



1/2 



(1) 



As only a small amount of radiation is transferred at these 
angles, we exclude them from our code and extrapolate the 
radiation field properties into this angle range. 

We calculate the opacities kj (cm 2 g -1 ) for polarization 
modes j=l (X-mode) and j = 2 (O-mode). The opacity kj is 
proportional to the imaginary part, Kj , of the complex refrac- 
tion coefficient: 



kj — 



Attv 
cp 



(2) 



where v is photon frequency and p plasma density. We present 
equations for Kj below. 

An example of the computed continuum opacities, together 
with the opacities at the cyclotron line and its harmonics, is 
shown in Figure [TJfor various 9. Notice that, at some angles, 
the opacities at the fundamental resonance and first two har- 
monics are larger than the continuum opacities in both modes. 
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Here hv — E is the photon energy, and v p is the electron 
plasma frequency, 

2 

v v = — ■ (5) 



The dimensionless energy shift x s relative to the central reso- 
nance energy in the s-th harmonic (s = 1 corresponds to the 
fundamental resonance) is defined by the following expres- 
sion 

_ £ - sb + s 2 p 2 b 2 /A x 

x s Tj 

x 4 



(6) 



where 



h hvB u j? t hv 



(7) 



= hv 



m e c 2 kT 



1/2 



|cos0|. 



In the same regime, the O-mode opacity is given by the 
following equation: 



K 2 : 



/3sin 4 6^ 



v 2 4(1 + cos 2 0)| cos 0| 
(3 + tan 2 - 2V) 2 
7r(l + cos 2 9) 2 (H 2 +F 2 ) 



b , b 
— tann — h 
4 2 



(8) 



cos 



exp(-x 2 ). 



Here H s = H(x s ) and F s = F(x s ) are the real and imagi- 
nary parts of the plasma dispersion function W(x s ): 



Figure 1. Opacities in cyclotron lines and continuum opacities in two polar- 
ization modes (top and bottom panels) in a magnetized plasma with T = 1 . 5 
MK, p = 1 g cm~ 3 , and B = 7 X 10 10 G, for four values of the angle 6 
between the magnetic field and the photon wave vector: 5° (thickest black 
curves), 30° (thick red curves), 60° (thin blue curves), and 80° (thinnest ma- 
genta curves). Continuum opacities for two values of 9 are shown by dotted 
curves (thick black curves for = 5° and thin magenta curves for 80°). 



This allows one to expect significant absorption lines at these 
resonances in the emergent spectra of magnetized NS model 
atmospheres. 

2.1. Opacities at the fundamental electron cyclotron 
resonance 

According to PSY80, the equations for the normal mode 
opacities near the fundamental resonance are qualitatively dif- 
ferent in the regimes of weak and strong vacuum polarization 
by magnetic field, i.e., at small and large values of the param- 
eter V(3 cos sin 2 0, where 



V 



3 x 10 28 cm 



B 



(3) 



n e is the electron number density, B the magnetic field 
strength, and Bq = 4.414 x 10 l3 G is the magnetic field 



at which E r 



m e c . In the case of CCOs, with their rela- 



tively weak magnetic fields, the regime of weak vacuum po- 
larization is applicable throughout the NS atmosphere. In this 
regime, the vacuum polarization does not affect the X-mode 
opacity: 



vl 1 



cos 2 



4/3 1 cos 0| 



\pK exp(— x\). 



(4) 



W(z) = H(z)+iF(z) = - 



■ du. 



(9) 



H(z) is the well-known Voigt function, proportional to the 
Doppler profile in our approximation, 



H(x + iy) = — 



7T J-oc, - u ) 2 + y 2 

while F(z) is its Hilbert transform: 



du — 5- e 



F(x + iy) = — 

IT 



e u (x — u) 
(x — u) 2 + y 2 



1 Z' + OO 

y->0 1 / e 



du 



■ du. 



(10) 



(11) 



2.2. Opacities in higher harmonics 

The values of Kj in the harmonics s > 1 are given by the 
following equation (signs "plus" and "minus" correspond to 
X-mode and O-mode, respectively): 



^A(X±Y)eM-xt), 



where 



X = l + cos 2 + Csiii 2 0, 
Y = (1 + G 2 y 1/2 [2| cos0| + G(l - C) sin 2 0)] 



G 



sin 2 v% + Vv 2 {v 2 - 
2|cos0| v%v 

C 2 /3 2 sin 2 b 
C= Asb tmh 2' 



(12) 

(13) 
(14) 

(15) 
(16) 
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3. ATMOSPHERE MODELS 

We calculate the magnetic hydrogen atmosphere models us- 
ing our recently developed code (Suleimanov et al. 2009; 
SPW10), assuming the magnetic field normal to the stellar 
surface. We also assume that the atmosphere is fully ion- 
ized because, in the relatively low magnetic fields considered 
here, and temperatures of a few million Kelvin, the fraction 
of bound species in a hydrogen atmosphere is very low (see, 
e.g., Pavlov et al. 1995). Moreover, all the hydrogen lines 
and photoionization jumps lie in far (or extreme) ultraviolet, 
unobservable because of the interstellar absorption (e.g., the 
energy of the strongest hydrogen line atB = 10 11 Gis about 
26 eV, while the photoionization energy of the ground state is 
about 78 eV; see Roesner et al. 1984). Therefore, atomic lines 
are unobservable in CCOs, and bound species do not make a 
strong effect on the temperature run and the emergent spec- 
trum. 

We calculated two sets of atmosphere models with the same 
parameters as in SPW10. The surface gravitational accel- 
eration g = 10 14 cm s~ 2 for all the models. In the first 
set we fixed the magnetic field strength (B = 7 x 10 10 G) 
and computed models with three different effective temper- 
atures, T c ff = 1, 1.5, and 3 MK. The emergent flux spec- 
tra and temperature structures for the models computed with 
and without the collisionless cyclotron opacities are shown 
in Figure [2] A common property of the emergent spectra 
for the models with the cyclotron opacities are prominent 
absorption cores at the cyclotron energy and its harmonics, 
E s = sE c , e = 0.81 s keV (s = 1,2,.. .), in addition to the 
wide absorption wings caused by the quantum oscillations in 
the magnetic free-free opacities in the cold plasma approx- 
imation (SPW10). The equivalent width w s of the s-th ab- 
sorption feature decreases with increasing s as well as with 
increasing T c h (i.e., with decreasing the effective quantiza- 
tion parameter, b c ff = E c ^ e /kT e s ss 9.4, 6.3, and 3.1 for the 
three models). However, these dependences are weaker than 
in the models computed in the cold plasma approximation. 
For instance, wi= 243 (240), 213 (200), and 213 (180) eV at 
T off = 1, 1.5, and 3 MK, respectively, while w 2 = 234 (190), 
205 (140) and 146 (55) eV at the same temperatures, where 
the numbers in parentheses are equivalent widths for the mod- 
els calculated in the cold plasma approximation. 

The second set consists of four models with different mag- 
netic field strengths (B = 1, 4, 7, and 10 xlO 10 G) at the 
same effective temperature, T c g = 1.5 MK. The compari- 
son of the emergent spectral fluxes for the models calculated 
with and without thermal cyclotron opacities is shown in the 
top panel of Figure [3] Again, the equivalent widths of the 
absorption features decrease with increasing harmonic num- 
ber and decreasing magnetic field strength. The difference 
of the emergent spectra computed with and without thermal 
cyclotron opacities is similar to that seen in Figure [2] In par- 
ticular, the new models show substantially deeper cores of the 
absorption features at the cyclotron energy and its harmonics, 
while the continuum spectra remain virtually the same. 

The bottom panel of Figure [4] shows spectra of emergent 
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Figure 2. Top panel: Emergent flux spectra for magnetic NS atmospheres 
(B = 7 X 10 10 G) in the cold plasma approximation (thin solid lines) and 
with the collisionless cyclotron opacities included (thick solid curves) for 
three effective temperatures, T c g = 1, 1.5, and 3 MK. For clarity, the spec- 
tra for T e ff = 3 and 1 MK are shifted along the ordinate axis by factors 
10 + and 10" 1 , respectively. The dotted curves show blackbody spec- 
tra for the same temperatures. Bottom panel: Temperature structures for the 
same models. Solid curves correspond to the models with collisionless cy- 
clotron opacities included, dashed curves to the models in the cold plasma 
approximation. Notice the overheating of the outer layers in the models with 
the collisionless cyclotron opacities included. 

specific intensity for the model with T G fr =1.5 MK and 
B = 7 x 10 10 G, for various angles 9 between the line of 
sight and the normal to the surface. The quantum oscillations 
together with deep cyclotron cores in the intensity spectra are 
seen at any angle. The cyclotron cores depend on 9 stronger 
than the quantum oscillation wings. They are virtually invis- 
ible at small angles, in accordance with the disappearance of 
thermal cyclotron resonances at 9 — > in the O-mode for 
s = 1 and both modes for s > 1 — see Equations (8) and 
(17), and they become narrower and deeper with increasing 
9. The equivalent width of the absorption feature at the fun- 
damental resonance depends on 9 non-monotonously. For in- 
stance, wi ~ 228 (220), 180 (160), 380 (320) and 400 (330) 
eV for 9 = 5°, 30°, 60° and 80°, respectively (the equiva- 
lent widths calculated in the cold plasma approximation are 
in parentheses). 

In the top panel of the same Figure [4] the angular distribu- 
tions of the emergent intensity at five photon energies (near 
the fundamental resonance, the first harmonic, and three con- 
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Figure 3. Top panel: Emergent spectral fluxes for atmospheres with T e fj = 
1.5 MK and magnetic fields B = 1, 4, 7 and 10 xlO 10 G, calculated with 
and without including the collisionless cyclotron opacities (thick and and thin 
solid lines, respectively). For clarity, the spectra for £?io = 10, 7 and 1 are 
shifted along the ordinate axis by factors 10+ 2 , 10+ 1 and 10 _1 , respec- 
tively. Bottom panel: Temperature structures for models with and without 
the cyclotron opacities (solid and dashed curves, respectively), for the same 
parameters. For clarity, the temperature structures for B±o = 10, 4 and 1 are 
shifted along the ordinate axis by factors 3, 0.5 and 0.5, respectively. 



tinuum energies) are shown. Unlike the absorption features 
arising in the cold plasma approximation (see SPW10), the 
specific intensities at the cyclotron cores are stronger peaked 
towards the surface normal (i.e., towards the magnetic field 
direction) than those in the continuum. 

Spectra of emergent specific intensity in the O- and X- 
modes, for the model with the same parameters as in Figure 
IU are shown in Figure [5] The intensity spectra in the X-mode 
depend on 9 only slightly (bottom panel), while the intensity 
spectra in the O-mode, which dominates at E > E c ^ e , vary 
with similar to the intensity summed over polarizations (cf. 
Figure 4). 

The bottom panels of Figures [2] and [3] show the tempera- 
ture structures for the six atmosphere models in comparison 
with the temperature structures of the corresponding models 
computed in the cold plasma approximation. The main effect 
of the thermal cyclotron transitions is the temperature rise in 
outer atmospheric layers, caused by an additional cyclotron 
heating due to the resonance in the O-mode opacity. 

The depths of the spectrum formation (the column densities 



Figure 4. Top panel: Angular distribution of the emergent specific inten- 
sity at the fundamental resonance, first harmonic (s = 2), and three con- 
tinuum photon energies (indicated in the plot) for a NS atmosphere with 
B = 7 X 10 10 G, T cff = 1.5 MK. Bottom panel: Spectra of the emer- 
gent specific intensity for the same model at different angles to the surface 
normal (indicated in the plot). 

which correspond to t % e — 1) for the new model are shown in 
Figure [6] We see that in the region of cyclotron harmonics the 
O-mode comes from deeper (hotter) layers than the X-mode. 
Since these layers are much deeper than the outer layers over- 
heated by the cyclotron absorption (see Figures 2 and 3), the 
overheating does not affect the emergent radiation summed 
over polarizations. 

We also computed the radiative acceleration g ra d, which is 
defined by the following equation 



2tt 



5rad 



2 

E 



dp 



{ki + (ji)nr v {n)dn, (is) 



where /i = cos 9. The calculated <? ra d values are smaller than 
the surface gravity for all the computed model atmospheres 
(Figure|7]i and can be ignored. 

The vacuum polarization by the magnetic field might be 
significant in outer atmosphere layers, where the electron 
number density is very low. We have checked, however, that 
including the vacuum polarization into the calculations only 
slightly changes the temperature structure and the radiative 
acceleration at the very surface (Figure [8]), and it does not 
change the emergent spectrum at all. In the top panel of 
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Figure 5. Spectra of the emergent specific intensity in the O-mode (top 
panel) and X-mode (bottom panel) for the same atmosphere model as in Fig- 
ure[4] at different angles to the surface normal (indicated in the plot). 

the same figure, the radiation pressure distributions in the at- 
mospheres are also shown, and the agreement with the equi- 
librium thermodynamic value in the optically thick layers is 
demonstrated. 

4. DISCUSSION 

We have presented the first computations of fully ionized 
hydrogen atmospheres of NSs for magnetic fields B ~ 10 10 - 
10 11 G without using the cold plasma approximation. For 
the magnetic fields chosen, the electron cyclotron energy, 
E c>e ~ 0.1-1 keV, is within the range of energies where the 
thermal emission from isolated NSs is usually observed. We 
have shown that the thermal effects lead to formation of deep 
Doppler cores in the absorption features at the cyclotron en- 
ergy and its harmonics, and to a sharper angular dependence 
of the emergent intensity at the cyclotron features. The main 
contribution to the equivalent widths of the cyclotron features 
still comes from the feature wings, which are due to the quan- 
tum oscillations in the magnetic free-free opacity and can be 
considered in the cold plasma approximation. However, the 
contribution of the Doppler cores is significant, and it in- 
creases with increasing the quantization parameter b e g and 
the harmonic number. The equivalent widths of the features 
reach ~ 150-250 eV in the examples considered; they grow 
with increasing b e s and are lower for higher harmonics. 

In our previous work (SPW10) we suggested that the har- 
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Figure 6. Photon energy versus flux formation depth (column density corre- 
sponding to t 1 e =1) in the O- and X-modes for the atmosphere with T e ff = 
1.5 MK, £?io = 7. In the top panel the energy band around the fundamental 
resonance is shown in more detail. 

monically spaced absorption features in the X-ray spectrum 
of the CCO IE 1207 (Sanwal et al. 2002, Bignami et al. 
2003) are due to the quantum oscillations in the magnetic 
free-free opacity. We ignored the relativistic effects, such 
as the thermal cyclotron processes, because the thermal and 
cyclotron energies are essentially nonrelativistic for IE 1207 
(kT e f{/m e c 2 < E c . e /m e c 2 ~ 10~ 3 ). However, the detailed 
study presented here has demonstrated that these effects are 
important and must be included in the model atmosphere cal- 
culations. The reason is the large ratios of the cyclotron res- 
onance opacities to the continuum opacity (e.g., ~ 10 4 and 
~ 10 1 for the s = 2 and s = 3 harmonics, respectively, for 
the atmosphere parameters relevant to the case of IE 1207). 
Therefore, the inclusion of the thermal cyclotron opacities 
in modeling atmospheres of IE 1207 and similar objects is 
needed. 

In future work we plan to compare the model spectra ob- 
tained by integrating the specific intensities over the visible 
NS surface (at different rotation phases and for various ori- 
entations of the rotation and magnetic axes) with the phase- 
dependent spectra of IE 1207 to infer the properties of this 
NS. The integration should lead to a broadening of cyclotron 
features because of the nonuniformity of the surface magnetic 
field, without dramatic changes of equivalent widths. We note 
that the pulsations detected in the radiation of IE 1207 are 
strong only in the spectral features, and we believe that the 
strong angular dependence of the modeled emergent intensity 
in the cores of cyclotron harmonics (see Figures 4 and 5) can 
explain this fact and provide a potentially powerful tool for 
the determination of the NS parameters. 
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Figure 7. Relative radiation accelerations in all the computed model atmo- 
spheres (see Figs.|2]and[3j. The atmosphere parameters are indicated in the 
panels. 

In SPW10 we discussed the importance of quantum oscil- 
lations in the magnetic free-free opacity for the interpretation 
of the X-ray spectra of some other CCOs (J 1852+0040 in the 
Kesteven 79 SNR and J0822-4300 in the Puppis A SNR) and 
low-field radio pulsars. It is clear that the thermal cyclotron 
processes may also be significant for these objects and must 
be included in calculations of their model atmospheres. 

In the X-ray spectra of XDINSs, whose magnetic field are 
2-3 orders of magnitude higher than in CCOs, proton (or ion) 
cyclotron lines are probably observed. The thermal cyclotron 
processes should be less important for ions because the pa- 
rameter £ p = max(fcT, E CiP ) / m p c 2 is at least three orders of 
magnitude smaller than £ e , and the Doppler width of the ion 
cyclotron resonance is w 40 times smaller. However, the ra- 
tio of the ion resonance cyclotron opacity to the continuum 
opacity is also large for atmospheres of these NSs (> 10 7 , 
see, e.g. Figs. 2 and 4 in Suleimanov et al. 2009). Therefore, 
a contribution of the ion cyclotron harmonics to the emergent 
XDINS model spectra cannot be excluded and have to be in- 
vestigated 7 . 

In the model atmospheres computed here, the radiation 
transfer calculations were performed using the LTE approx- 

7 Here we mean the thermal Doppler cores of the ion cyclotron harmonics, 
not the small quantum peaks in the magnetic free-free opacities (see Potekhin 
2010 and references therein). 



Figure 8. Comparison of temperature structures (bottom panel) and relative 
accelerations due to radiation pressure (top panel) for models with (dashed 
curves) and without (solid curves) vacuum polarization taken into account. 
In the top panel, the corresponding radiation pressures are also shown. 

imation. Therefore, the source functions at the cyclotron line 
and its harmonics were taken equal to the Planck function, and 
the photon scattering in the cyclotron resonances was ignored. 
For a more accurate consideration, one should include non- 
coherent photon scattering in the cyclotron resonances. This 
is not simple (see, for example, Nagel 1981), and it has never 
been considered self-consistently. We plan to develop magne- 
tized atmosphere models with non-coherent photon scattering 
in the cyclotron line and its harmonics in future work. 
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